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ABSTRACT 

This paper presents a sample of "cold front" clusters selected from the Chandra archive. The clus- 
ters are selected based purely on the existence of surface brightness edges in their Chandra images 
which are modeled as density jumps. A combination of the derived density and temperature jumps 
across the fronts is used to select nine robust examples of cold front clusters: 1ES0657-558, Abell 1201, 
Abell 1758N, MS1455.0+2232, Abell 2069, Abell 2142, Abell 2163, RXJ1720. 1+2638, and Abell 3667. 
This sample is the subject of an ongoing study aimed at relating cold fronts to cluster merger activity, 
and understanding how the merging environment affects the cluster constituents. Here, temperature 
maps are presented along with the Chandra X-ray images. A dichotomy is found in the sample in 
that there exists a subsample of cold front clusters which are clearly mergers based on their X-ray 
morphologies, and a second subsample which harbor cold fronts, but have surprisingly relaxed X-ray 
morphologies, and minimal evidence for merger activity at other wavelengths. For this second subsam- 
ple, the existence of a cold front provides the sole evidence for merger activity at X-ray wavelengths. 
We discuss how cold fronts can provide additional information which may be used to constrain merger 
histories, and also the possibility of using cold fronts to distinguish major and minor mergers. 
Subject headings: galaxies: clusters: individual: 1ES0657-558, Abell 1758N, MS1455.0-h2232, 

AbeU 2069, Abell 2142, Abell 2163, RXJ1720. 1-^2638, AbeU 3667, Abell 665, 

Abell 2034 — X-rays: galaxies: clusters 



1. INTRODUCTION 



The hierarchical nature of structure formation is 
well understood th e oretically jP rcss & Schcchtcr 1974; 
ILacev fc Coli Il99l IPeeblesl[T993l: ISoringel e t al. 200ffl 
and well supported observationaliy (iSnioot et al. 19921: 
Hernauist et all 119961 : iCole et all 120051 : iJeltema et all 



2005f ). At the present epoch, clusters of galaxies are the 



largest and most massive large-scale structures that are 
virialized and they are growing hierarchically through the 
gradual accretion of surrounding matter. Occasionally, 
two clusters of roughly equal mass fall together in a major 
merger, releasing around 10^'^ erg of gr avitational poten- 
tial e nergy during the merger process (jMarkevitch et all 
[T999h . To a de gree, the majority of clusters exhibit evi- 
dence for merger activity. However it is the scale of the 
cluster mergers, i.e., the distinction between major merg- 
ers and the continuing minor mergers/infall of smaller 
subsystems, that is difficult to quantify. A gauge of 
merger activity capable of distinguishing major and mi- 
nor mergers would greatly assist to understand mergers 
and their effects on clusters and their constituents. 

A solution may be provided by high resolution 
X-ray observations w ith the Chandra observatory 
(jWeisskopf et al.l l2002f) . Among its first discoveries 
was the detection and characterization of the "cold 
front" phenomenon in the intraclust er media (ICM) 
of Abell 2142 (Markcv itch et al.l [200l and Abell 3667 
(jVikhlinin et al.l l2001at ). Cold fronts reveal themselves 
as edge-like features in the Chandra images and are 
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well modeled as contact discontinuities which form at 
the interface of cool, dense gas and hotter, more dif- 
fuse gas. Simulations show cold fronts can arise in 
more than one way during c luster mergers (jPoole et al.l 
120061 : lAscasibar fc Markevit ch 2006). Broadly speak- 
ing, cold fronts can be separated into two general 
types: "remnant core" and " sloshing" (for a review see 
iMarkevitch fc Vikhlininl[2007f) . 

Remnant core type cold fronts were first invoked to 
explain the origin of the cold fronts in Abell 2142 
(jMarkevitch et al]l2000f ). In this case, the discontinuity 
occurs between the interface of an infalling subcluster's 
cool core and the hotter ambient ICM of the main clus- 
ter. The subcluster's core is bared to the hotter ICM 
after ram pressure stripping removes its less dense, lower 
pressure outer atmosphere during its passage through 
the main cluster. The cool core survives longest be- 
cause, at first, it remains bound to the infalling dark 
matter core and, after separating from that, because 
of its high density. This explanation for Abell 2142's 
cold fronts has since been abandoned in favor of 
a "sl oshing" type mechanism (jMarkevitch fc VikhlininI 
l2007f l. An excellent example of remnant core type 
cold front is observed in the Bullet cluster (1ES0657- 
558; IMarkevitch et all 120021 ) while a similar mecha- 
nism acts to form cold fronts in the elliptical galax- 
ies NGC 1404 and NGC 4552 (where the cold fronts 
are at the inte rface of the galaxy's gas and the cluster 
ICM; Machac ck et all 12005. 2006). Furthermore, rem- 
nant co re type cold fronts have also bee n seen in both ide- 
ahzed (lAscasibar fc Markevitch . 20061: iPoole et al l 12006 



iSpringel fc Farraii 120071: iMastropietro fc Burkertl l2008f) 



and cosmologic al ("Bialek et al. 20021: 'N agai fc Kravtsovl 
l200llOnuoraet al. 2003i: lMathis~et al. 2005f l' simulations 
of cluster mergers. 
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Gas "sloshing" came to the fore as an explanation for 
the existenc e of a cold front in th e "relaxed" cluster 
Abell 1795 (|Markevitch et al.ll200fl ). These cold fronts 
occur when the stably stratified ICM is disturbed in such 
a way that the cool, dense gas residing in the cluster 
core is displaced from its position at the bottom of the 
gravitational potential well. The cold front is the con- 
tact discontinuity formed at the interface of the low en- 
tropy core gas and the higher entropy gas at larger radii. 
In most postulated scenarios, the disturbance is gravita- 
tional m nature with its origin in an infalling subclustcr 
Ce.g.lTittlev fc HenriksenilllOStrAscasibar fc Markevitd^ 



[2006ir However, there have been suggestions the distur- 
bance may be a hydrodynamic ph enomenon in the form 
of weak shocks or acoustic waves fChurazov et al.| |2003l : 
IFuiita et al1l2004D . while [Markevitch ct al. (2003b) sug- 
gest AGN explosions may be responsible, although the 
example of Hydra A used was later shown to be bet- 
ter inter preted as a weak shock (iNulsen et al.ll2005D . To 
date, the I Ascasibar fc MarkevitchI ( 20061 ) studv provides 
the most comprehensive set of simulations for the slosh- 
ing scenario and show it is possible that the cold fronts 
observed in otherwise relaxed appearing clusters (based 
on their X-ray morphologies) can be produced by slosh- 
ing induced in relatively minor mergers. Currently, there 
is no observational evidence that the cold fronts in these 
relaxed appearing clusters are related to merger activity. 

While cold fronts offer interesting insights int o the 
physical properties of the IGM (e.g. , JVikhlinin ct al. 
2001b; .Churaz ov fc I nogamovl [200l [L vutikov 2006; 
Markevitch fc V ikhlininI 120071 : IxTang et a l. 2007), it is 



how their presence relates to cluster merger activity, 
along with the tantalizing possibility of utilizing cold 
fronts as a gauge for cluster mergers, which provides the 
motivation for the study presented in this paper. The 
utilization of cold fronts as merger gauges becomes even 
more compelling when considering those clusters which 
appear to harbor relaxed X-ray morphologies — the exis- 
tence of a cold front provides the only indication of ongo- 
ing merger activity at X-ray wavelengths. Much work has 
gone into studying cold front formation with simulations, 
but what has been lacking to date is a systematic study 
of the relationship between cold fronts and other dynam- 
ical indicators of cluster merger activity. The number of 
clusters observed by Chandra to harbor cold fronts has 
grown significantly and now allows a representative sam- 
ple to be assembled in order to perform such a study. To 
this end, this paper focuses on the selection of a sample 
of clusters with cold fronts in their Chandra X-ray im- 
ages. The selection criteria are designed to ensure that 
the high quality X-ray data shows an unambiguous cold 
front in each member of the sample. 

To date, our goal of establishing whether there is a 
link between the presence of cold fronts in clusters and 
evidence of recent dynamical growth has been met, at 
least partially, in two c ases, Abell 1201 (lOwers et al.l 
[20b9b) and Abell 3667 (jOwers et al.] I2009al ). In these 
studies, significant dynamical substructure was detected 
using spatial and redshift information provided by large 
samples of spectroscopically confirmed cluster member 
galaxies. In these clusters, the presence of a cold front 
is directly related to the presence of substructures, and 
thus ongoing merger activity. Future work will endeavor 
to go one step further and study the galaxy properties 



within the cold front clusters. 

The outline of this paper is as follows. In Section[2]the 
procedure for selecting cold front clusters is presented. In 
Section [3] the technique used to generate projected tem- 
perature maps for each cold front cluster is described. 
In Section [4] the Chandra X-ray images and temperature 
maps for the clusters in the cold front sample are pre- 
sented, along with a brief description of the X-ray prop- 
erties and evidence for merger activity from the litera- 
ture. In Section m the results and sample are discussed. 
In Section [5] summary and conclusions are presented. 

2. SAMPLE SELECTION CRITERIA 
2.1. Initial Selection 

To ensure selection of only the best defined cold front 
clusters, the following criteria were set. The data must 
be of high enough quality that statistically robust mea- 
surements of density and temperature can be made on 
either side of the edge in surface brightness for charac- 
terization. The bulk of the cluster X-ray emission must 
be captured on the Chandra chips, whilst resolution and 
cosmological surface brightness dimming affects must be 
minimized. Clearly, the cluster must also exhibit an edge 
in the surface brightness, indicating the existence of a 
density discontinuity caused by a cold front. 

The initial sample was selected from those clusters with 
publicly available data within the Chandra archive (as of 
July 2006) and the preceding conditions were imposed as 
the following quantitative criteria: 

• a total Chandra ACIS-I and/or ACIS-S exposure 
time exceeding 40 ks; 

• cluster redshift in the range 0.05 < z < 0.3; 

• a significant edge in the X-ray surface brightness. 

The density discontinuities caused by cold fronts are 
generally quite visible in the Chandra images, and their 
existence is easily verified by the observation of edges 
and compression of isophotal contours overlaid onto the 
images. Thus, for each Chandra data set meeting the 
first three selection criteria, the 0.5 - 7keV band im- 
age was scrutinized for edge-like features and regions of 
compressed isophotes. It is noted that the morphology 
of the X-ray emission was not considered during the se- 
lection process, nor were any clusters selected based on 
prior knowledge of a cold front in their ICM. The sig- 
nificance of the density discontinuity was quantified by 
modeling the surface brightness across the edge with a 
broken power law density model (see Appendix E| ■ The 
surface brightness edge must also be reasonably sharp, 
i.e., in order to be classified as a front, they must be 
reasonably well fitted by the density jump model. 

This model was used to further cull the sample and is 
presented in Appendix [Xj A description of the method 
used to fit the model to the data follows. 

2.2. Fitting the surface brightness across the edges 
Corresponding to the broken power law density model. 



the surface brightness model given in Equation IA4I was 
incorporated into the Sherpa fitting software, which is 
part of the Chandra Interactive Analysis of Observations 
[CI AO; version 3.4 is used here) package. Point sources 
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were identified using the CIAO WAVDETECT tool and 
removed for fitting. The images were restricted to the 
energy range 0.5 - 7keV, so that background and instru- 
ment effects were minimized, whilst still allowing enough 
source photons to obtain good fits to the models. Where 
possible, separate observations had their data and ex- 
posure maps co-added prior to fitting. The background 
count rates for Chandra observations is known to change 
with time, and also to differ between the ACIS-S and 
ACIS-I arrays. For these reasons, observations taken 
at significantly different times, or taken using different 
chips were not co-added, but were simultaneously fit- 
ted with separate background components. The fit was 
restricted to regions where there are edges (the blue re- 
gions shown in the top left panels of Figures [2] through 
fTU]) . and the initial inputs for the center of the spheroid, 
the ellipticity, e, the orientation angle, 9, and the radius 
at which the edge occurs, i?/, were obtained by over- 
laying an elliptical region that best approximated the 
surface brightness contours across the edge. The cen- 
ter of curvature was not necessarily the position of the 
X-ray centroid of the cluster, and the semimajor axis of 
the ellipse, 9, was oriented to bisect the edge. Since the 
center, 9 and Rf are highly degenerate, the center and 
9 were fixed to those values determined from the man- 
ually fitted ellipse and only Rf was determined directly 
from the fit. A fixed, uniform background component 
was added to the model during fitting. The background 
was determined by fitting the entire Chandra observed 
area using one or two beta models plus a constant com- 
ponent. The beta model/s account for the cluster emis- 
sion, and the constant component for the background. 
We note that this constant background component does 
not correctly account for the unvignettcd instrumental 
background. This does not affect our results because: a) 
the edges are found in regions where the cluster emission 
is dominant and b) the quantity of interest is the break 
in the surface brightness profile which is insensitive to 
both the background level and to small changes in the 
background slope (when the background does not domi- 
nate). Each model was multiplied by the exposure map, 
generated using the CM MKINSTMAP, ASPHIST and 
MKEXPMAP tools, which account for instrumental ef- 
fects such as vignetting, Quantum Efficiency (QE), QE 
nonuniformity, bad pixels, dithering and effective area. 

After fitting the density model (Appendix IA4p . the 
sample was limited to include only those clusters where 
the density jump exceeded a factor of 1.5 at the lower 
limit of the 90% confidence interval. The confidence in- 
tervals were estimated using the PROJECTION tool in 
Sherpa, which varies the parameter of interest along a 
grid of values whilst allowing adjustable model parame- 
ters to settle to their new best-fitted values. This crite- 
rion eliminated several clusters with clearly visible sur- 
face brightness edges from the sample. Thus, any de- 
tectable cold fronts in clusters meeting the other selec- 
tion criteria above that were missed in the search of the 
archive are unlikely to have density jumps exceeding a 
factor of 1.5. The 11 clusters which meet the above crite- 
ria, along with the parameters for the fits, are presented 
in Table [TJ The corresponding surface brightness distri- 
butions with the surface brightness profiles for the best 
fitting density models overlaid are plotted in Figure [TJ 
Interestingly, Abell 2142 is the only one of the 11 clusters 



to require a significant ellipticity in the fits. It is noted 
that density jumps determined from the surface bright- 
ness fits do not depend on the curvature of the front (e^; 
in Appendix I A4p . but they do depend on the assumption 
that the density profiles can be expressed in terms of the 
same coordinate (r in equation I Al[) on both sides of the 
front. The Chandra Observation Identification number, 
cluster center and redshift are presented in Table [2] along 
with global spectral properties (Section 12.3. ip . To con- 
firm that these structures are cold fronts, temperature 
measurements on each side of the front are required (see 
Section [1X2D. 

2.3. Preparation of data for Spectral measurements 

Since the ACIS chips on board Chandra record both 
spatial and energy information for incident photons, it 
is possible to extract spatially resolved spectra from the 
images. As for all astrophysical spectral observations, 
the data need to be processed carefully in order to ex- 
clude events not related to the source (e.g. cosmic rays), 
calibrate the energy and also to account for the back- 
grounds. Here, the procedure used to prepare the data 
for spectral extraction is described. 

The archived Chandra pipeline data were reprocessed 
starting with the level 1 event files and using CIAO ver- 
sion 3.4 with the Calibration Database version 3.4.0. Ob- 
servation specific bad pixel files were produced and ap- 
plied, as were the latest gain files. Observations taken 
at focal plane temperatures of — 120°C had corrections 
applied for charge transfer inefficiency and time depen- 
dent gain. Where observations were taken in VFAINT 
telemetry mode, the more thorough VFAINT mode of 
cleaning was used for improved rejection of cosmic rays. 
Only events with ASCA grades of 0, 2, 3, 4 and 6 were 
retained during the analysis. 

Backgrounds were taken from the blank sky observa- 
tions appropriate for the epoch of observation*. Since the 
blank sky backgrounds are taken from observations with 
low Galactic foregrounds and low soft X-ray brightness, 
the soft X-ray flux was checked in the vicinity of each 
cluster using the ROSAT all sky R45^ count rates, and 
compared to the blank sky background rates. If the R45 
fluxes were anomalous (e.g. MS1455.G+2232, AbeU 2163 
and Abell 3667), a region relatively free from cluster 
emission was selected and spectra extracted to model 
the soft component using a s ingle or double MEKA L 
model at zero redshift (see e.g. iMarkevitch et al.ll2003al ). 
The soft component was included in all spectral analyses, 
with the normalization weighted by a geometric factor to 
allow for the different extraction areas. The background 
and source data were processed using the same calibra- 
tion files, bad pixel files and background filtering, with 
the backgrounds being reprojected to match the obser- 
vations. 

To match the blank sky backgrounds, the observations 
were filtered for periods of high cosmic ray background 
caused by flares. The point sources and regions contain- 
ing cluster emission were excluded and a light curve was 
extracted and binned in time intervals of 259 seconds in 
the energy range 0.3-10 keV for the ACIS-I array, whilst 
the ACIS-S array data were binned in intervals of 1037 

See http://cxc. harvard. edu/contrib/maxim/acisbg/j 
^ See ^http: / /heasarc. gsfc.nasa.gov/docs / tools.html ^ 
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Fig. 1. — Surface brightness profiles across the edges (black points with error barsVwith the surface brightness for the best fitting density 
model plotted in red. The corresponding model parameters can be found in TablelT] The clusters shown are (from left to right, top row 
through bottom row) fES0657-558, Abell 665, Abell f201, Abell f758, MSf455.0+2232, Abell 2034, AbeU 2069, AbeU 2f42, Abell 2f63, 
RXJ1720. 1+2638 and Abell 3667. 



seconds in the energy range 2.5-6 keV. The hght curve 
was analyzed in Sherpa using the lc_clean.sl script and 
only periods where the background count rate was within 
20% of the quiescent rate were included. 

During the spectral analysis of Abell 1758N, it was 
necessary to account for a soft flare present for the du- 
ration of the obse r vation . The method described in 
iDavid fc KempneH (|2007 ) was used, whereby the flare 
shape was modeled in PHA space in a region relatively 
free from cluster emission, using the cutofFpl model in 
XSPEC with power law index = 0.15 and an exponential 
cut off at E=5.6 keV, with the normalization allowed to 



vary. The flare model was then incorporated into spec- 
tral fits with all parameters fixed and the normalization 
set to its best fitted value from above multiplied by a 
geometric factor correcting for the different extraction 
region areas. 

2.3.1. Global Temperatures 

To obtain a picture of the overall thermal properties 
of the 11 clusters, global temperatures and metallicities 
were measured. Where possible, a region covering the 
majority of the cluster emission was selected and spec- 
tra were extracted using the CIAO DMEXTRACT tool. 
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TABLE 1 

Parameters of the density model fits for the 11 clusters with significant density discontinuities. 



Cluster 


(x.y) 


e(A6») 




e 


^1 






02 




(J2000) 


(deg) 


(kpc) 




(10-8)=' 









1ES0657-558 


104.58518, -55.941469 


(310 - 50) 


•^•^-0.2 





22.7+^-^ 


0+°-^ 


2, 


g+0.4 
■"-OA 




Abell 665 


127.74663, 65.842256 


252 (233 - 270) 


192+1 





''■"-O.l 


^■'-0.4 


0, 


O-I-0.3 
■^-0.3 


1 o+°-2 

^•'-'-0.2 


Abell 1201 


168.2197, 13.4509 


254 (240-267) 


5031* 





2 0+"-3 


J^-U-0.3 


0, 


4+0.2 
■^-0.1 


1 q+0.3 
^^•■^-0.2 


Abell 1758N 


203.1603, 50.5575 


98 (45-150) 


150+2 





3 n+0.4 


-0 i+°-i 


1, 


1+0.2 

■-^-0.2 


1 1+°1 


MS1455. 0+2232 


224.3125, 22.3427 


132 (114-150) 


150^1 





13.2+i* 


1 9+0. 2 
^■^-0.2 


3, 


4+0.9 


1 6+°-3 


Abell 2034 


227.53396, 33.495151 


120 (104-126) 


521+2 





„ 7+0.2 
'-0.2 




0. 


7+0.1 
■'-0.1 


9 a+0-2 


Abell 2069 


231.12201, 29.994898 


25 (3 - 46) 


3831* 







1 q+0.5 
^■■'->0.2 


0, 


o+O.l 

■■-■-o.i 


1 S+0-2 
^•^-0.2 


Abell 2142 


239.5863, 27.226989 


40 (10 - 100) 


360ti 






q+°-i 


1, 


O+O.l 

■''-O.l 


-I O+0.04 


Abell 2163 


243.97508, -6.1117719 


315 (310-320) 


856^2 





18.3«-i 


-0 4+°-^ 


5, 


9+0.4 
■^-0.4 


o o+O.l 

■^•«-o.i 


RXJ1720. 1+2638 


260.03982, 26.624822 


234 (205 - 263) 


190+1 





10 i+°-8 


4: + °-^ 


2, 


4+0.3 


1 4+0.1 


Abell 3667 


303.14776, -56.819424 


230 (222 - 238) 


4^^+0.3 





4 4+0.2 


-0 7+°-'^ 
^- ' -0.1 


0, 


fj+0.1 

■"-0.1 


1 5+°-2 



Note. — For a description of the density model parameters Ai, A2, 
^ Units of Ai and A2 arc photons/ cm'^ / s / arcsec^ 

with point sources detected with WAVDETECT excised. 
The regions within which the spectra were extracted are 
shown in red in the X-ray images presented in the top 
left panels of Figures [2] through [lOl The CIAO tool MK- 
WARF was used to create auxiliary response files (ARE), 
which account for spatial variations in quantum efficiency 
(QE), effective area and also the temporal dependence of 
the effective area due to molecular contamination on the 
optical blocking filter. For observations performed with a 
focal plane temperature of -120° C, the CMC MKACIS- 
RMF tool was used to create redistribution matrix files 
(RMF), while MKRMF was used for this purpose for 
observations performed at — 110°C. Since the responses 
depend on the chip position, they are calculated in re- 
gions in chip coordinates and are weighted by the 0.5-2 
keV brightness of the cluster within that region, with the 
final ARE and RMEs being count weighted. For clusters 
with multiple observations where a more recent, signifi- 
cantly longer exposure was available, the earlier, shorter 
observations were excluded from the global temperature 
measurement. 

The spectra were fitted in the 0.5-9.8 keV range for 
ACIS-S observations and 0.6-9.8 keV for ACIS-I obser- 
yations in XSPEC j Arnaud 1996) us ing a MEKAL model 
(iMewe et alj 119851 ,1986^ .KaastTal Il992t iLiedahl et all 
119951 ). for a hot, diffuse, single temperature plasma, mul- 
tiplied by WABS, a photo-electric absorption model ac- 
counting for Galactic absorption, where the neutral hy- 
drogen column de nsity is initially fixed to t he Galactic 
value provided bv lDickev fc LockmanI (|1990[ ). The spec- 
tra were binned so that each energy bin contains at least 
30 counts and the normalization, temperature and abun- 
dance were fitted in the first iteration, with the column 
density fixed to the Galactic value. The metal abun- 
dances a re measured relative to th e solar photospheric 
values of lAnders &: Grevessd (|1989[ ). Where the data al- 
low, the column density was then freed and the data 
re-fitted, with the best fitting results along with their 
associated 90% errors presented in Table [2j For the clus- 
ters 1ES0657-558, AbeU 2142, Abell 2163 and Abell 3667 
the measured iVn values are inconsistent with the Galac- 
tic values within the 90% confidence limits, although 
they are consistent with values measured by others in 



oi and Q2, sec Appendix IA] 

the literature (e.g. iGovoni et al.|[200l : iMarkevitch et al.l 
120001: iVikhlinin et al.ll2001aD . In the clusters Abell 665 
and RXJ1720. 1+2638, the TVh values are also inconsis- 
tent with Galactic values, however the temperature and 
abundance measurements appear to be relatively un- 
affected and they are consistent w i th previously mea- 
sured values (Maz zotta et al.ll200ll IGovoni et all 12004 
for RXJ1720. 1+2638 and Abell 665, respectively). The 
LUMIN function in XSPEC was used to determine the 
integrated, unabsorbed X-ray luminosity from the best 
fitting model in the energy range 0.5-9.8 keV, with the 
results presented in Table [2| 

2.3.2. Temperatures across the edges 

Cold fronts are expected to have very different thermal 
properties when compared to shock fronts. The pressure 
{P ~ TiekT) across a cold front should be continuous, 
whilst a shock front should show a pressure discontinu- 
ity. The specific entropy, specified here in terms of the 

entropy index, E = kTrie should also change rapidly 
across a cold front since, regardless of the physical mech- 
anism causing the cold front, the front occurs at the in- 
terface of cool, dense low entropy gas and the hotter, am- 
bient higher entropy ICM. Thus, a jump in the entropy 
across a front is expected which will also be opposite to 
the entropy change expected across a shock front. To 
ensure the observed surface brightness edges are indeed 
cold fronts, the temperature was measured just inside 
and outside the edges and used in combination with the 
density jumps measured above to derive changes in pres- 
sure and entropy across the edges. 

The precision of the pressure measurements across the 
fronts is determined by limitations in both the data and 
the assumptions used in modeling the surface brightness 
across the edges to obtain densities. For example, in ob- 
taining densities it has been assumed that the gas density 
is a function of the elliptical radius. This is an approx- 
imation, which may be poor if, e.g., there is emission 
from gas projected onto the edge region, or there are 
rapid, nonradial variations in the gas distribution. These 
effects are reduced by confining attention to the most 
prominent fronts, where the X-ray emission is generally 
dominated by gas close to the front. Furthermore, while 
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TABLE 2 

Clusters with density jumps exceeding 1.5 at the lower 90% confidence interval. The exposure times presented in 

PARENTHESES ARE THE CLEANED EXPOSURE TIMES USED TO GENERATE THE TEMPERATURE MAPS. THE A^h VALUES PRESENTED IN 

PARENTHESES ARE THE GALACTIC VALUES. 



Cluster 


Obsid 


Exposure 


RA 


DEC 


Redshift 


kT 


Z (Z0) 


ix 








{Kb) 




{O ^KJKJKJ ) 








^iu erg/s; 


xlO^o 


1ES0657-558 


3184, 4984, 4985, 


551 (535) 


06:58:27 


-55:56:47 


0.296 


14 7+0-3 
-0.3 


r, 09 + 0.03 


43.0 


3.7t;5-j(6.53) 




4986, 5355, 5356 














5357, 5358, 5361 


















Abell bbo 


531, 3586 


40 (35) 


r\o on A cr 

08:30:45 


65:52:55 


0.182 


8-5-0.7 


n qn+0.10 


13.6 


"7-1-1. 2/^ 

2-7_o.8(4-24) 


Abell 1201 


4216 


40 (22) 


11:13:01 


13:25:42 


0.168 


c; 0+0.3 




4.00 


1.61 


Abell 1758N 


2213 


58 (42) 


13:32:32 


50:30:36 


0.279 


^•^-0.5 


0.56i»;i? 


14.2 


1.05 


MS1455.0+2232 


534, 4192 


102 (90) 


14:57:15 


22:20:30 


0.258 


4 5+°-i 
^•^-0.1 


'-'■^'-'-0.04 


14.6 


3.1 


Abell 2034 


2204 


54 (53) 


15:10:13 


33:31:42 


0.113 


0-^-0.2 


0.26tHi 


5.4 


1.58 


Abell 2069 


4965 


55 (39) 


15:23:58 


29:53:26 


0.116 




0.28t»;0^ 


4.6 


1.96 


Abell 2142 


1196, 1228, 5005 


68 (45) 


15:58:16 


27:13:29 


0.089 




41+0 02 
'^■^J^-0.02 


17.5 


3.3i°;^(4.20) 


Abell 2163 


545, 1653 


81 (80) 


16:15:34 


-06:07:26 


0.201 


1.: 0+0.8 




34.1 


lbAt°ol (12.1) 


RXJ1720.1+2638 


304, 549, 1453 


60 (50) 


17:20:09 


26:37:35 


0.164 


6 

"•J^-o.i 


41+0 05 


14.0 


3.2t0;«(4.06) 




3224, 4361 












Abell 3667 


513, 889, 5751 


534 (473) 


20:12:34 


-56:50:26 


0.055 


7-1+0.1 


'^■^'J-o.oi 


3.9 


i.22+_li (4.71) 



5752, 5753, 6292 
6295, 6296 



pressure continuity holds in the immediate vicinity of 
the cold front, if there is a stagnation region at the front 
then in upstream regions away from the cold front, where 
the velocit y gradient is high, the pressure will vary (see 
Figure 6 of IVikhlinin eralll2001ar) . Thus, pressure mea- 
surements must be made in regions as close to the front 
as possible. The size of these regions is most affected 
by the temperature measurements, which require a large 
number of counts in order to be accurate, meaning the 
precision of the pressure measurement across the edges is 
limited mostly by the temperature measurements, rather 
than the shortcomings of the surface brightness model. 

In order to minimize the above effects, regions from 
which spectra were extracted for temperature measure- 
ments are chosen to lie as close as possible to the fronts. 
The opening angles for these regions correspond to those 
used for measuring the surface brightness profiles in Sec- 
tion 12.21 (see the blue regions plotted in the left panel 
of Figures [2] through [T0|) . Three exceptions exist in 
Abell 1201, AbeU 665 and MS1455. 0+2232 where for the 
region just outside the edge, a slightly larger region with 
a larger angular extent was used (MS1455. 0-1-2232 has a 
slightly larger radius compared to the region where the 
surface brightness profile was measured, too), so that 
there arc enough counts to measure the temperature 
with sufficient accuracy. Using a larger angular extent 
was preferred to increasing the radial extent for the rea- 
sons discussed above. Regions were chosen with a mini- 
mum of 700 counts in the 0.5 - 7keV energy range after 
background subtraction. The green arcs shown in top 
left panels of Figures [2] through [10] show the inner and 
outer limits for the regions used to measure the tempera- 
tures on the inside and outside of the edges, respectively. 
Where a green arc is not plotted, the radial limit used 
corresponds to that used in measuring the surface bright- 
ness profile (i.e. the blue regions). 

The procedure for fitting the temperature inside and 
outside the edge was as follows. First, the temperature in 
the region outside the edge was fitted using an absorbed 



MEKAL model with the column density and metallicity 
fixed^ to the values derived in Section B.3.11 Second, the 
spectrum from the region just inside the edge was fitted 
with two absorbed MEKAL models. The first MEKAL 
component accounts for the denser parcel of gas lying 
within the edge. The second MEKAL component ac- 
counts for gas lying in projection along the line of sight 
which is assumed to have the same thermal properties 
as the gas fitted in the region outside the edge. Thus, 
the temperature of this second MEKAL component was 
fixed to that measured for the region outside the edge. 
In XSPEC, the normahzation of the MEKAL model is 
defined as 

where Da is the angular diameter distance to the cluster, 
rig and Up are the electron and proton densities. This is 
proportional to the emission measure, so the normaliza- 
tion of the second MEKAL component was fixed to the 
value derived for the outer region, however it was cor- 
rected by a factor accounting for the different emission 
measures expected from the different volumes probed. 
Using the density model for (r) in Equation lAll for 
Rp < r < Rout (where Rout is the outer radius of the re- 
gion used to measure the temperature outside the front, 
shown in each cluster's Chandra image [left panel of Fig- 
ures [2] through [10]), along with the parameters presented 
in Table]!] Equation]!] can be integrated for both regions 
inside and outside the edge and the ratio used to obtain 
the correction factor applied to the normalization of the 
second MEKAL component. 

Table ]3] lists the results of the temperature measure- 
ments across the edges, as well as the density, entropy 

^ We have tested the effect of allowing the metallicity to vary 
during the fitting procedure using the Abell 3667 observations. We 
find that, within the errors, the temperature measurements are 
not affected, despite measuring metallicities of Zi = 0.84^q'22 and 
Z2 = 0.49^027 inside and outside the front, respectively. 
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TABLE 3 

Properties across the edges in surface brightness. 



Cluster 


Density Jump 


kT Inside 


kT Outside 


Pressure Jump Entropy Jump kT2 for shock 






(feTi) 







1ES0657-558 
Abell 665 
Abell 1201 

Abell 1758N 
MS1455.0+2232 
Abell 2034 
Abell 2069 
Abell 2142 
Abell 2163 
RXJ1720. 1+2638 
Abell 3667 



3.0(2.7-3.4) 
2.2(1.6-3.2) 
2.1(1.7-2.7) 
1.9(1.7-2.2) 
2.0(1.7-2.4) 
2.3(2.1-2.5) 
2.0(1.7-2.4) 
2.0(1.9-2.1) 
1.9(1.8-2.0) 
2.1(1.9-2.3) 
2.6(2.4-2.8) 



5.8 (5.3-6.3) 
9.2 (7.1-13.1) 

3.6 (2.9-4.7) 

7.9 (5.7-12.1) 

3.7 (3.3-4.3) 
6.6 (4.9-10.2) 
4.5 (3.5-6.3) 

7.5 (6.1-9.5) 

8.6 (7.1-10.8) 
5.4 (4.4-6.8) 

3.7 (3.5-4.0) 



19.7(16.9-22.2) 

8.6 (6.6-11.9) 

5.7 (4.0-9.4) 
10.3 (7.8-15.4) 

7.0 (5.6-9.2) 

4.8 (3.9-6.3) 
6.0 (4.4-9.1) 

17.6 (13.1-25.4) 
13.2 (9.9-18.9) 
7.4 (6.3-9.3) 
8.2 (7.3-9.3) 



0.9 (0.6-1.3) 

2.3 (1.0-6.3) 

1.4 (0.5-3.2) 

1.5 (0.6-3.4) 

1.0 (0.6-1.8) 

3.1 (1.6-6.6) 
1.5 (0.6-3.4) 
0.8 (0.5-1.5) 

1.2 (0.7-2.2) 
1.5(0.9-2.5) 
1.2 (0.9-1.6) 



0.1 (0.1-0.2) 
0.6 (0.3-1.4) 

0.4(0.2-0.8) 
0.5 (0.2-1.1) 
0.3 (0.2-0.6) 

0.8(0.4-1.6) 
0.5 (0.2-1.0) 
0.3 (0.2-0.5) 
0.4 (0.24-0.8) 
0.5 (0.3-0.7) 
0.2 (0.2-0.3) 



1.6 (0.9-2.4) 

4.8 (1.6-9.1) 

1.9 (1.0-3.2 ) 
4.8 (2.9-8.3 ) 

2.1 (1.5-3.0 ) 

3.2 (2.0-5.6) 
2.6 (1.6-4.4) 

4.3 (3.3-5.7) 
5.2 (4.1-7.0) 
1.8 (1.6-2.1) 

1.4 (1.2-1.8) 



and pressure jumps for each cluster in the sample. For 
comparison, the temperature expected on the less dense 
side of the edge is computed under the assumption that 
the density jump is caused by a shock front using the 
Rankine-Hugoniot condition 

m ^ (7 + i)-grt(7-i) 

kT, (7 + l)-^(7-l) ^' 

where 7 — 5/3 is the adiabatic index for a monatomic 
gas, quantities with subscript 1 refer to those on the in- 
side of the edge (i.e. denser side), and those with sub- 
script 2 refer to quantities on the outside of the edge 
(less dense side). Abell 2034 and Abell 665 both have 
best fitting kT2 values which are lower than kTi and, 
although the error bars would allow the temperature 
to be continuous across the front, the defining criterion 
of pressure continuity across a cold front is violated in 
these clusters (although for Abell 665 the lower error 
bar marginally allows for pressure continuity). A more 
likely interpretation for these edges in surface brightness 
is that of shock fronts, and this can be also seen by com- 
paring the expected temperatures for the pre-shock gas 
given by the Rankine-Hugoniot shock conditions with the 
temperature measured on the outside of the edge. Thus, 
Abell 2034 and Abell 665 are excluded from the cold front 
sample. However, shock fronts strongly imply the exis- 
tence of merger activity within these clusters and con- 
tamination in a cold front sample by shock fronts does 
not have a detrimental effect if the primary objective is 
to use the fronts as signposts for merger activity. 

The remaining nine clusters all have thermal properties 
across their surface brightness edges which are consistent 
with that expected for a cold front. The cold front sam- 
ple is comprised of these nine clusters and their X-ray 
images and projected temperature maps are presented 
in Section m 

3. TEMPERATURE MAPS 

Temperature maps can offer further insights into the 
overall dynamical state of a cluster, revealing regions 
which have been heated by shocks, dissipation of tur- 
bulence and adiabatic compression during a merger, 
and are also useful for comparing to the temperature 
structures found in simulations of cluster mergers. For 
the cold front clusters, temperature maps were cre- 
ated using the broad energy band method described in 



iMarkevitch et all (j200Clf ). Briefly, for each cluster, source 
images were produced in several energy bands in the 
range 0.5 - lOkeV. The energy bands were chosen to 
include at least 5000 counts per band (or 10000 where 
the data allow). Exposure maps correcting for mirror vi- 
gnetting, QE variations, photoelectric absorption by con- 
taminant on the filters and exposure time were produced 
for each energy band. Background images in the corre- 
sponding energy bands were produced using the blank 
sky observations described in Section 12.31 and were nor- 
malized by the ratio of the source to background 9-12 
keV counts and subtracted from the source images which 
were then divided by the exposure maps. Each image 
was smoothed using the same variable width Gaussian, 
where the width of the Gaussian, cr(r), is smallest in the 
brightest regions, i.e. the regions with the most source 
counts, and becomes larger as the brightness lowers. Dif- 
ferent initial values of cr(r) were trialled until a value was 
found which allows statistically significant temperatures 
to be measured (i.e., the la errors did not exceed 30% 
of the best fitting value at each pixel) whilst maintaining 
a good degree of spatial resolution. The noise in each 
pixel was determined from the raw, uncorrected images 
and weighted accordingly to allow for the effects of the 
smoothing. 

The temperature maps were produced by fitting an 
absor bed MEKAL model to the band fluxes for each 
pixel (|Markevitch et ani2000f ). The absorption column 
was set to either the Galactic value, or the best fitting 
value from Section [2.3.11 The metal abundance was set 
to the best fitting average cluster value derived in Sec- 
tion [23111 The model was multiplied by an on- axis ARF 
to correct for the energy dependent on-axis mirror effec- 
tive area, including the chip quantum efficiency (QE). 
A flux-weighted spectral response matrix was generated 
from a large cluster region, and was binned to match the 
chosen energy bands. 

Where there were multiple observations for a single 
cluster, images and exposure maps were co-added to pro- 
duce single images in each energy band. However, if 
there was a significant time gap between observations, 
or the observations were taken on different chips, sepa- 
rate images were produced for each observation in each 
energy band and the data were simultaneously fitted to 
produce a single temperature map. Regions in the tem- 
perature map where the la errors were greater than 30% 
of the best fitting temperature, or where the best fit- 
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ting temperature exceeds 25 keV, were excluded from the 
maps which are shown in the top right panel of Figures [2] 
through [TOl Regions near the edges of the maps, where 
the counts are low, are also excluded. Furthermore, for 
Abell 1758N, the spatial variations of the soft flare are 
not accounted for in the generation of the temperature 
map. However, the temperature structure seen in the 
right panel of Figure [6] agrees qualitatively with structure 
seen in maps which have been corrected for this spatial 
nonuniformity^. As we are interested in qualitative tem- 
perature structure, we chose not to include modeling of 
the flare spatial structure here for simplicity and note 
that our conclusions are not affected. 

4. NOTES ON INDIVIDUAL CLUSTERS 

Here, the X-ray images, temperature maps and r-band 
optical images for each of the cold front clusters are 
presented and discussed. There is a clear dichotomy 
in the sample when considering X-ray morphology in- 
somuch that the sample can be divided into those clus- 
ters which appear disturbed (i.e. they are clear merg- 
ers with a complex appearance and/or substructure) and 
those which, aside from the existence of a cold front, ap- 
pear relaxed. To highlight this, the sample is partitioned 
based on X-ray morphology with Section 14.11 contain- 
ing the disturbed clusters and Section containing the 
relaxed appearing clusters. Within these sections, the 
clusters are arranged in decreasing order of density jump 
strength. For each cluster, we discuss further evidence 
for merger activity present in the literature and, where 
possible, relate this to the existence of the cold front. 

It is noted that all of the clusters in the sample have 
previously published work based on Chandra data (aside 
from Abell 2069), and some of these studies included 
temperature maps of the cold front clusters. However, 
for the purpose of comparison, it is best to present tem- 
perature maps which have been produced in a homoge- 
neous manner. In addition to this, many of the clusters 
have since been reobserved by Chandra (for example, 
Abell 2142 and Abell 3667) meaning the maps produced 
here are of a higher statistical quality and resolution than 
those previously published. Where pertinent, the tem- 
perature maps presented here are compared to the pre- 
viously published versions. For each temperature map, 
the color scale used was chosen so that different colored 
regions reflect regions where temperatures differ by at 
least ~ ItT. 

4.1. Clusters with Disturbed X-ray Morphologies 
4.1.1. 1ES0657-558 - Bullet cluster 

1ES0657-558 was discovered serendipitously as an 
extended X-ray s ource in the Einstein Slew Survey 
() Elvis et al.l I1992D and was later found to be one of 
the h ottest, most luminous clusters known ([Tucker et al.l 
Il998| ) with redshift z = 0.296 and a veloc i ty dis - 
persion of Uy ~ IJOOkms-ijB arrena "eFall [200l) . 

iMarkevitch et all ()2002[ ) presented the first Chandra ob- 
servations of 1ES0657-558, which revealed a high average 
temperature of ~ 15keV and a clearly disturbed mor- 
phology consisting of a main cluster, along with a sub- 
cluster to the west and two sharp surface brightness edges 

^ See http: / /cxc. harvard. edu 7ccr/proceedings/03-proc/prese-| 
ntations / markevitch2 / index, html 



on the western side of the subcluster. These features are 
seen in greater detail in the Chandra image presented 
in the top left panel of Figure [2] which makes use of a 
much deeper exposure. The inner edge is a cold front 
and the pressure across it is continuous within the er- 
rors, whilst the outer edge is a shock front with a large 
pressur e jump consistent with a shock with Mach number 
M = 3 (IMarkevitch &: VikhHni"nll2007t IMarkevitch et all 
120021) . This is one of the best examples of a remnant 
type cold front, where the cool core of the merging sub- 
cluster survives the initial infall, forming a cold front at 
its interface with the hot, shocked main ICM. The less 
dense outer subcluster ICM has been stripped from the 
denser core which survives because of its high density. 
This stripped gas can be seen in the X-ray image as a 
fan-like structure extending out behind the subcluster. 

The top right panel of Figure shows the tempera- 
ture map with surface brightness contours overlaid. The 
temperature structure is complex showing pockets of 
hot (~ 20keV) gas mixed with cooler (~ 12keV) gas. 
The only regular structure lies in the region containing 
the cold and shock fronts, where it can be seen that 
the subcluster/cold front coincides with the region of 
^ 7keV gas (slightly hotter than the temperature mea- 
sured above of ~ 5.8 keV due to projection effects) , signif- 
icantly cooler than the surrounding gas, whilst the shock 
front coincides with an arc of 18 — 20keV gas. The hot 
regions on the eastern side of the main cluster are prob- 
ably shock heated regions of either stripped subcluster 
gas, or main cluster ga s . Previous temper ature maps 
(IMarkevitch et al.l 120021 : iGovoni et all \2004) also show 
1ES0657-558 to have a complex morphology with shock 
heated gas of up to 20 keV, consistent with the map pre- 
sented here. However, the temperature map we present 
here has higher resolution and statistical quality, since it 
is derived from almost an order o f magnitude more Chan- 
dra data. iMillion fc AllenI (|2008f ) used the same data set 
to produce a temperature map using an adaptive binning 
technique. The temperature map presented here agrees 
well with the Million fc Allen temperature map. 

Th e clust er hosts a luminous radio halo (Liang et al.l 
12000(1 ■ and IGovoni et all (|2004[ ) showed that the halo 
brightness peaks correspond to hot spots in their tem- 
perature map, and the halo extends along the shock 
front. 1ES0657-558 has been the subject of a thorough 
campaign to map its mass distribution using the com- 
bination of weak a nd strong lensing (Clowe ct al. 20061; 
iBradac et al.l [20061 ) . the results of which have shown an 
offset between the lensing mass distributions and the X- 
ray centroids providing compelling evidenc e for the ex- 
istence of dark matter ([Clowe et al.l I2006D. These re- 
sults, comb ined with the observatio ns of iBarrena et al.l 
dlOOZ) and IMarkevitch et all ([2002f ) show that the sub- 
cluster is probably caught just after its first (slightly off 
center) core passage, traveling from the south-east, and 
the merger occurs mostly in the plane of the sky, with 
the velocity difference betwee n the main cluster a nd the 
subcluster only ^ 600kms~^('Barrcna et al."2002') . The 
veloci ty of the shock derived by Markevitch & Vi khlininI 
([2007f ) is 4700 km s~^, although in reality the subcluster 
is unlikely t o be traveling at such a high velocity (see, 
for ex ample, iSpringel fc Farraill2007tlMilosavlievic et al.l 
[2001 . 
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Fig. 2. — Top left: background subtracted, expos ure cor rected, 0.5 - 7keV Chandra image of 1ES0657-558. The red region shows the area 
used to measure the global temperature in Section 12.3.11 and the blue region is where the surface brightness profile was fitted. The green 
arcs show the radial limits used for measuring the temperature inside and outside the cold front (the cold front delineates the two regions). 
North is towards the top of the page, and East is to the left. Top right: this image shows the projected temperature map produced as 
described in Section |3] The color scale is tuned so that different colors correspond to regions where temperatures differ by approximately 
1(T. Bottom left: Magellan IMACS r-band image (courtesy of D. Clowe). Bottom right: zoomed in version of the temperature map showing 
the region containing the cold front. The top right, bottom left and bottom right images have X-ray surface brightness contours overlaid. 
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4.1.2. Abell 3667 

Abell 3667 was one of the original cold front 
clusters observed with Chandra and, along with 
Abell 2142, provided the d e fining characteristics of 
cold fronts (Vikhlinin_gt_alJ 1200 lal) . Abell 3667 
has an X-ray luminosity of Lx(0.4 — 2.4 keV) = 
5.1 X lO'^'^e rgs"^ (lEbeling et al.' '19961 ), temperature 



kT ~ 7keV (iKnopp et al.lll9 96: Markev itcheFaLl [l998l : 
IVikhlinin et al.ll2001at IBriel et al.ll2004l ). Abell richness 
class 2 and a redshift of z=0.055. In the top left panel of 
Figure [3] the Chandra X-ray image of the central core re- 
gion of Abell 3667 is presented revealing a complex mor- 
phology which is elongated along a south-east to north- 
west direction with the cold front being the most promi- 
nent feature to the south-east. The front appears regular 
at small opening angles, but fans out at larger angles to 
form a mushroom-cap-like shape. The same morphol- 
ogy is also prominent in the temperature map in the 
top right panel of Figure [3] which shows clearly that the 
front delineates regions of cool (~ 4.5 keV) gas from the 
ambient ICM wit h temperature (~ 7keV). As noted by 
IBriel et all (I2004f). the g a s at th e tip of the front is coolest 
4keV). E emz eFaLl (pool ) use simulations to show 
motions inside a remnant gas core are induced due to 
the flow of gas around the cloud. These internal motions 
move lower entropy gas towards the tip of the cold front 
and, as the gas moves to the tip of the cold front, it cools 
due to adiabatic expansion, enhancing the temperature 
contrast at the tip of the front. The temperature map 
presented here is in good agreement with those presented 
in the hterature (iMarke vitch et" al]|1999 : IVikhlinin et al.l 
I2001at iMazzotta et al.l l2002: Briei.eLalJ[2Q0^, although 
the longer combined exposure time of the Chandra data 
presented here gives a more statistically significant map 
with higher spatial resolution. The map shows a com- 
plex morphology of hot and cool gas patches, including 
a patch of significantly hotter (~ 9 — lOkeV) gas to the 
north-west of the cluster center. 

At the cold front, a density jump of 2.6 ± 0.2 is mea- 
sured (see Table [3]). This value is incons istent with th e 
value of 3.9 ± 0. 8 measured bv.Vikhlinin et all (|2001aD . 
This is because IVikhlinin et al.l 's density measurement 
just outside the front is obtained by extrapolating their 
/3-model electron density model from larger radii to the 
front position. Their aim was to determine the pres- 
sure in the undisturbed region of the flow around the 
front, so they tried to correct for gas within the "com- 
pression" region just outside the front where they note 
that the surface brightness ~ 70 kpc exterior to the front 
exceeds that predicted by their /3-model. For the density 
measurement made here, the power law density model 
describing the surface brightness outside the front is ap- 
proximately limited to this ~ 70 kpc region, meaning a 
higher density j ust outside the front is measured com- 
pared to that of IVikhlinin et all , which leads to a lower 
density jump. 

Abell 3667 presents clear evidence for merger activity 
at both optical and radio wavelengths. Observations at 
radio wavelengths show Abell 3667 is one of a handful 
of objects which harbo r two diffuse radio relic sources 
(|Rottgering et al.lll997f ) — one to the north-west and one 
to the south-east. The radio observations also show two 
tailed radio galaxies, one of which lies on an axis joining 



the two domin ant cluster galaxies wit h its tail pointing 
along this axis (|Rottgering et al.lll997l ). Radio relics are 
thought to trace merger shock fronts and lRoettiger et all 
(1999) were able to reproduce the basic radio morphology 
of the relics using magnetohydrodynamic/N-body simu- 
lations of a slightly off-center merger with a 5-to-l mass 
ratio and the subcluster impacting from the south-east 
to the north-west, being observed ~ 1 Gyr after core pas- 
sage. The simulations also reproduced the basic observed 
X-ray morphology and galaxy distribution. 

Optical observations have shown the galaxy den- 
sity distribution is bimodal, with the second den- 
sity enhanceme nt centered on th e second br i ghtest 
cluster galaxv (iProust et al.l 119881 : iSodre et al.l 119921 : 
lOwers et al.ll2009a|) . Signiflcant emission associated with 
the secondary density enhanc ement was de t ected in the 
ROSAT X-ray observations of'Knopp et all (11996'), lying 
just out of the field o f view of the Chand ra observations. 
Using weak lensing, iJoffre et"al] (|2000[ ) found evidence 
for substructure in their mass maps whi ch coincides 
with this secondary density enhancement. ISodre et al.l 
(jl992| ) found the second D galaxy has a velocity offset of 
~ 400 kms~^ compared the the central D galaxy, and the 
velocity distribution surrounding this second D galaxy 
appears to indicate that it forms a dynamical subunit 
within the cluster. 

Abell 3667 was the subject of comprehensive optical 
follow-up with the aim of showin g the c old front was 
related to merger activity in lOwers et al.l (|2009af) . Us- 
ing combined spatial and velocity information for a large 
sample of sp ectroscopically confirmed cluster members, 
lOwers et al.l found the cluster can be separated into two 
major compone nts, coinciding with the density enhance- 
ments found bv ISodre et al.l (|1992 ^ and separated in pe- 
culiar velocity by ^ 500kms~^, plus a number of smaller 
subgroups. These observations indicate Abell 3667 is un- 
dergoing a major merger which has produced the ob- 
served cold front. 

4.1.3. Abell 1201 

Abell 1201 is one of the least well studied clust ers in 
the s ample and has an A bell richness class 2 ( Abell et al.l 
IT989I ). redshift z=0.168 (jStruble fc Rood 11999,) with an 
X-ray luminosity Lx (0 -1 - 2.4 keV) = 2.4 x lO^'^ergss-i 
(|Bohringer et alJ 120001). The top left panel of Figured] 
shows the Chandra image of Abell 1201 which is elon- 
gated along a south-east to north-west axis, with the 
main cold front clearly visible to the south-east, a bright 
central core which also has a cold front towards the 
north-west and a remnant core also visible to the north- 
west. The temperature map shown in the top right panel 
of Figure [4] reveals a hot ~ 7 keV region lying between 
the main cluster core and the remnant core, indicating 
some form of heating here probably caused by adiabatic 
compression due to core motion. There is no evidence for 
significant temperature differences at the location of the 
remnant core (c.f. 1ES0657-558 where the remnant core 
clearly stands out in the temperature map). This may 
be attributed to the heavy smoothing required to obtain 
well-defined temperatures, since this observation was sig- 
nificantly affected by strong flares, so that nearly half of 
the exposure was rejected during cleaning. There is a 
cool (~ 4keV) stream of gas connecting the cold front 
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Fig. 3. — Same as Figure[2]but for Abell 3667 and with a Supercosmos Sky Survey rp-band image at the lower left. 



to the cool core of the main cluster in the temperature 
map. This indicates that the cold gas associated with the 
front is gas which has been displaced from the cool core 
of the cluster, and the front is probably associated with a 
spiral-like structure (for example as seen in the temper- 
ature maps for MS1455.0-t-2232 and RXJ1720. 1+2638 in 
the top right panels of Figures [9] and [H respectively) seen 
in projection, similar to those seen in the x projection 
image presented in Figure 19 of lAscasibar fc MarkevitchI 
(l200l . 

Detailed follow up multi-object optical spectroscopy 
of Abell 1201 has been undertaken and is presented in 



lOwers et al.l ()2009bD . This study shows a bimodal galaxy 
density distribution with the main component centered 
on the main X-ray core and a second component coinci- 
dent with the remnant core to the north-west. The com- 
bination of spatial and velocity information also reveals 
the presence of substructure coincident with the remnant 
core and offset in peculiar velocity by ^ 430kms~^, in- 
dicating a merger is occurring with the majority of its 
velocity in the plane of the sky. The perturbation from 
this merger has given rise to the cold fronts caused by 
sl oshing of th e ICM 
lEdge et "all (|2003D discovered a small scale gravita- 
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Fig. 4. — Same as Figure|2]but for Abell 1201 and with an SDSS r-band image at the lower left. 

tional arc in Abell 1201, and proposed the observed Abell 2069 is a richness class 2 (|Abell et aLlll989| ) clus- 
arc properties could be explained by a mass distribu- ter at redshift z=0.116, which was originally shown to 
tion which was either highly elliptical or bimodal. The have multiple condensations in the gas and galaxy distri- 
data allowed a bimodal model with a secondary mass butions by Gio ia et al.l (|1982f ) and is part of a large super - 
clum p lying to th e south-east of the cluster center, how- cluster containing some 6 clusters (jEinasto et al.l[T997[ ). 
ever. lEdge et al.l fav ored the model with high ellipticity. The top left panel of Figure [5] shows the Chandra im- 
Given the results of lOwers et al.l ()2009bl ). it would seem age of Abell 2069 (binned to 7.8" pixels) which shows 
that a bimodal mass distribution is the most likely ex- the main compon ent (Abell 2069A in the terminology of 
planation for the data, although the center of mass of iGioia et al.lll982! ) to be elongated along a south-east to 
the second component is coincident with the north-west north-west axis. There are two bright elliptical galaxies 
remnant core. aligned along this axis which differ in p eculiar velocity 

by only ~ 230 km s"^ (|Gioia et al.l[T98g ). Abell 2069B 

4.1.4. Abell 2069 can be seen to the north-east and it hosts a cold front 
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towards the north-west along with a fan like extension of 
X-ray emission towards the south-east. There is a domi- 
nant elliptical galaxy associated with Abell 2069B which 
has a redshift z=0.1178 and is surrounded by an over- 
density in the p rojected galaxy number surface density 
(jGioia et al.l 119821 . Assuming the redshift of the domi- 
nant elliptical is representative of Abell 2069B, then its 
peculiar velocity with respect to the mean cluster red- 
shift is only 485 kms~^. Although the cold front is asso- 
ciated with the Abell 2069B, it is difhcult to determine 
whether the nature of this cold front is consistent with 
a contact discontinuity occurring at the interface of the 
lower entropy Abell 2069B ICM and the higher entropy 
Abell 2069A ICM, or a cold front induced by sloshing of 
the Abell 2069B gas. Given the small peculiar velocity 
offset of Abell 2069B and Abell 2069A, it can be con- 
cluded that the majority of the subcluster motion is in 
the plane of the sky. 

The resolution of the temperature map presented in the 
top right panel of Figure[5]is limited by a lack of photons. 
However, it does show that the projected temperature 
around Abell 2069B is significantly cooler than the ICM 
surrounding it, and also that Abell 2069 A hosts a com- 
plex temperature distribution with regions of 6 — 8 keV 
gas interspersed with pockets of cooler ^ 4keV gas. The 
temperature structure in Abell 2069A may be due to a 
core crossing by Abell 2069B. In particular the region 
of 6 — 8keV gas to the east may be due to the north- 
east bound passage of Abell 2069B. In this case, the cold 
front in Abell 2069B is more likely to be due to slosh- 
ing of the Abell 2069B gas. However, the existence of 
the two dominant elliptical galaxies may be indicative of 
merger activity within Abell 2069A itself and this could 
also explain the complex temperature structure. 

4.1.5. Abell 1758N 

Abell 1758 is a rich (Abell richness class 3; 
lAbell et al.l Il989f) dou ble cluster at redshift z=0.279 
(|Struble fc RooTfl999l) . The Chandra observations. 



which were originally presented in lDavid_& Kempner 
were taken on the ACIS-S3 chip and contain only 
the northern component of the double cluster which is re - 
ferred to as Abell 1758N. The lDavid fc Kempneii (|200i ) 
study showed the X-ray morphology is irregular with two 
main components placed on an axis running from the 
south-east to the north-west. The iDavid &: Kempneii 
■|2[)( ) 1 ;■ study also showed that the morphology of the 
north-west component is fan-like with a tail of gas ex- 
tending to the south, and a cold front located towards 
the north, suggesting the motion of the substructure is 
towards the north. These features can be seen in top 
left panel of Figure [6] where Chandra X-ray image is 
reproduced. Using a truncated sph eroid to model the 
surfac e brightness across the front, IDavid &: Kempneii 
(|2004l ) measured a density jump of 1.6 ± 0.2 (la er- 
rors) which is consistent within the errors with the value 
presented in Table [3l The north-west and south-east 
structures are coincident with peaks in t he projected 
galaxy number density maps presented in iDahle et al.l 
and also w i th sub structure in the lensing maps of 
lOkabe fc Ui^Ietsul (|2008l ). There is a second less signifi- 
cant front further south-east of the south-east structure 
for which IDavid fc Kempne^ (|2004[) measured a density 
jump of ~ 1.5 ± 0.2. They used a hardness ratio map 



to show that the gas inside the front is cooler than out- 
side, thus concluding it is also a cold front caused by the 
interface of low entropy subcluster gas with the higher 
entropy ambient ICM. Coincident with the parcel of low 
entropy gas responsible for the south-east cold front is 
another overdensity of galaxies housing a bright ellipti- 
cal galaxy. 

The top right panel of Figure [6| shows the temperature 
map for Abell 1758N, which only covers a small area 
due to the low number of photons detected, particularly 
in the outskirts. The map shows the two components 
host significantly cooler gas compared to the hotter gas 
in the outskirts. The absence of a sharp temperature 
jump at the cold front is mainly due to the poor effective 
resolution of the temperature map. This means the cold 
front regions are contaminated by the hotter gas from 
the outskirts, thus the temperature measured here is the 
average of the two contributions. There is a bridge of cool 
gas joining the two subclusters, and a tongue of cool gas 
trailing the south-east subcluster which is presumably 
ram pressure stripped gas due to the collision of these two 
subclusters. The features seen in the temperature map 
are qualitati vely similar to tho s e see n in the hardness 
ratio map of IDavid fc Kempned (|2004|). 



Bas ed on the X-ray characteristics. IDavid fc Kempnej 
(|2004f ) presented a merger scenario for Abell 1758N 
whereby the system is being viewed after a recent off- 
axis pericentric passage of two ^ 7keV clusters, and the 
significantly hotter gas in the outskirts has been heated 
by merger shocks. Indirect evidence for merger activ- 
ity occurs at radio wavelengths where, at the position of 
the south-east cold front, there is a narrow-angle tailed 
(NAT) radio galaxy which ha s its tail extending fro m the 
north-west to the south-east (jO'dea fc Owenl[l985f ). and 
its head coincident with a fainter gala xy to the nor t h-east 
of the dominant south-east elliptical. iBliton et al.l (|1998l ) 
showed that NAT galaxies are preferentially found in 
clusters with substructure, concluding that at least part 
of the relative velocities required to produce a NAT is 
due to bulk gas motions generated during a merger. Fur- 
ther Chandra observations with longer exposures would 
enable the measurement of pressure differentials at the 
cold front, thus constraining gas velocities and allowing 
the components of the relative velocities of the gas and 
NAT galaxy to be constrained. 

Given the fan-like morphology and its coincidence with 
a galaxy overdensity, it appears the north-west cold front 
in Abell 1758N is a merger induced cold front generated 
at the interface of colliding subclusters. 

4.1.6. Abell 2163 

Abell 2163 i s one of the hottest (kT = 15.5 keV; 
iMaughan eFall I2008D . most l uminous (LxjO.l - 
2.4keV] = 22.9 x lO'^^ gj.g g-i; lEbeling et all [199^ 
clusters of galaxies known. It has a redshift of 

z = 0.201, Abell richness class 2 (jAbell et al.l 

|1989| ) and has been the subje ct of many stud- 
ies at multiple wavelength s (e.g.; iM aurogordato et alj 
20081: [Feretti et al. 2004; Markeyiteh fc Vikhlinin 2Jigi[ 
Feretti et a"Ll |2001; Squires ct al. 1997; Holzapfcl et al.l 



1997f). The lMarkevitch fc Vikhhnin, (,2001.) Chandra ob- 
servations of Abell 2163 revealed a cluster with an ir- 
regular morphology with an extension along a direction 
going north-east to south-west, as seen in deeper Chan- 
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Fig. 5.— Same as Figure[2]but for Abell 2069 and with an SDSS 



dra image presented in top left panel of Figure [T] There 
is a cold front in the south-west which appears to have 
a mildly convex shape and is not as sharp as others seen 
in this sample (e.g. the cold front seen in Abell 3667), 
possibly indicating the front has significant inclination 
with respect to the plane of the sky or is in the process 
of breaking up. There is a group due north and a hint of 
excess emission east of the core emission. The tempera- 
ture map shown in the top right panel of Figure [7] shows 
the cluster exhibits a complex multi-temperature gas dis- 
tribution with the cold front associated with ^ 8 keV gas 
and surrounded by hotter ~ llkeV gas, consistent with 
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the temperature jump given in Table [3] within the er- 
rors. The regions corresponding to the X-ray extension 
towards the south-east and the excess X-ray emission 
to the east are significantly hotter than the rest of the 
cluster at > 14keV. The infalling group to the north 
is significantly cooler than the main cluster at 5-6 keV 
according to the temperature map. The temperature 
map presente d here is in good ag reement with the maps 
publis hed by iGovoni et al.l ((2004!) and Markevitch et al] 
(j2001h . although we obtain higher resolution and statis- 
tical quality by simultaneously fitting the two Chandra 
data sets to produce the temperature map. 
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As mentioned above, Abell 2163 has been studied 
extensively at muhiple wavelengths and evidence for 
merger activity presents itself in the majority of these 
observations. The first indications that Abell 2163 was 
an atypical cluster came from combined Ginga and Ein- 
stein X-ray obser vations which showed an extremely 
high temperature (|Arnaud et alj 1992). Further ROSAT 
and ASCA observations confirmed the high tempera- 
ture and also showed evidence for temperature and sur- 
face brightness struc ture in the ICM, suggesting evidence 
for merger activity (|Markevitch et all 119941 : lElbaz et al.l 
[1995; Mar kevitchlll996D . H igh resolution Ch a ndra tem- 
perature maps presented in lMarkevitch et all (|2001[ ) and 



iGovoni et all (|2004f ) revealed this temperature struc- 
ture in greater detail, showing in particular that the 
central region harbors complex temperature structure 
which is consistent with a major merger, along with the 
group towards the north seen in the ROSAT images. 
I Squires et all (|1997f) shed light on the dynamical state 
of the cluster showing the distributions of both mass 
and galaxies are bimodal in the center, with the peaks 
in each distribution roughly coincident, and concluded 
the cluster is not in a relaxed dynamical state. At ra- 
dio wavelengths, Abell 2163 hosts a luminous radio halo, 
a possible radio relic, along with 3 tailed radio galaxies 
all with tails oriented in the same direction towards the 
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Fig. 7. — Same as Figure|2]but for Abell 2163 and with a CFHT r-band MegaCam image at the lower left. 



west (iFferetti et al.ll2001[ ). Finally. iMauro gordat o et all 
(|2008f ) present an analysis of the cluster dynamics from 
a sample of 361 cluster member redshifts which show 
a multi-modal velocity distribution, strong velocity gra- 
dients along a north-east to south-west axis along with 
multiple clumps in the galaxy surface density. These re- 
sults were confirmed by Owcrs ( 2008) where the analysis 
of a sample of 491 cluster member redshifts showed a 
significantly skewed velocity distribution and significant 
dynamical substructure when spatial and velocity infor- 
mation was combined. In particular, substructure was 
found to coincide with the cold front seen in Figure [71 
strongly supporting the argument that this cold front is 



caused by an ongoing cluster merger. 

4.2. Clusters with relaxed Appearing X-ray 
Morphologies 

4.2.1. RXJ1720.1+2638 

RXJ1720.1-)-2638 was the first of the otherwise re- 
laxed appearing clusters which were revealed to have cold 
fronts by Chandra (Mazzotta et al. 2001). It has an X- 
ray luminosity of Lx(0 .1 - 2.4 keV) = 7.3 x lO^'^ erg s~^ 
(iBohringer et al.ll2000f) and redshift z = 0.1605 (|Owersl 
I2008f) . The top left panel of Figure [H] shows the Chan- 
dra X-ray image of RXJ1720. 1-1-2638 which reveals a re- 
laxed overall morphology with a bright core and a cold 
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front towards the south-east. iMazzotta fc Giacintuccil 
(|2008) also report a cold front towards the north-west 
of the core which has a density .iump of ^ 1.68. It is 
noted that IMazzotta fc Giacintuccil ()2008f ) measured a 
density jump of 1.66 ± 0.05 for the south-east cold front, 
incon sistent with the value given in Table [S This is be- 
cause |Miiiott^^iGiacintucc3 ()2008D measured the den- 
sity jump over a larger opening angle, which included 
regions where the front is less well defined. The temper- 
ature map presented in the top right panel of Figure [8] 
reveals more about the nature of the cold front. There 
is a finger of cool gas starting at the surface brightness 
peak and spiraling in an anticlockwise direction before 
terminating at the cold front, also seen in the tempera- 
ture maps of lMazzotta fc Giacintuccil (|2008D , and similar 
to the temperature structure s een in the simulations of 
lAscasib ar fc MarkevitchI (|2006l . see their Figure 7). On 
larger scales the temperature map is patchy with regions 
o f hot ^ 8 keV gas s urrounded by ambient ^ 6keV gas. 

iDahle et al.l (|2002i ) report on the distribution of galax- 
ies, light and dark matter within RXJ1720. 1-1-2638, find- 
ing the light distribution to be fairly regular and dom- 
inated by the central cluster galaxy. The weak lensing 
mass distribution appears to be less regular, with the 
main signal roughly coincident with the cluster center 
and a secondary peak detected 2 — Sarcmin to the 
south-west. The galaxy number density also appears to 
be fairly regular, although there is an extension of the 
contours towards the south-west in the direction of the 
secondary mass clump detected in the weak lensing maps. 
RXJ1720.H-2638 is by all accounts a fairly relaxed look- 
ing system, and, along with MS1455. 0-1-2232, is thus a 
very important test case within our sample, given the hy- 
pothesis that cold fronts are a result of major merger ac- 
tivity. For this reason, RXJ1720. 1+2638 was selected for 
follow up observations using multi-object spectr o scopy . 
The results of the study are presented in lOwersI ()2008f ) 
and show evidence for substructure in the galaxy density 
distribution and also when spatial and velocity informa- 
tion is combined. However, the substructure's relation 
to the cold front remains ambiguous, and further obser- 
vations are required. 

4.2.2. MSI455. 0+2232 

MS1455.0+2232 (also known as ZW7146) was discov- 
ered serendipitously as an extended X-ray source in the 
Einstein Ob serv atory Extended Medium-Sensitivity Sur- 
vey (EMSS; iGioia et al.lll99df) and lat er reported to hos t 
a massive ISOOMgyr-i cooling flow (|Allen et al.lll996D . 
The top left panel of Figure M shows the Chandra X- 
ray image of MS1455. 0-1-2232 which has a relaxed X- 
ray morphology, harboring a slight ellipticity with ma- 
jor axis running from the south-west to north-east with 
a mild asymmetry along this axis. Also along this axis 
is the cold front as outlined by the blue annular seg- 
ment in the top left panel of Figure [H Inspection of 
the temperature map presented in the top right panel of 
Figure |9| reveals a spiral-like feature of cool gas begin- 
ning at the cool core and ending at the position of the 
cold fr ont. This feature is also s een in the temperature 
map of lMazzotta fc Giacintuccil (|2008[ ). which is derived 
from the same Chandra data set, and shows good over- 
all agreement with the map presented here. The cold 
front appears similar to those seen in the simulations of 



lAscasibar fc MarkevitchI (|2006f) (see Figure 7 and Figure 
19) where an off-center merger of a gasless dark matter 
halo disturbs the cool cluster core, offsetting both the 
gas and dark matter from its initial position. The offset 
low entropy gas does not fall back into the core along 
a purely radial path, since it has acquired angular mo- 
mentum from the infalling dark matter halo, and thus a 
spiral pattern is formed. 

Further evidence for merger activity in 
MS1 455. 0-1-223 2 com es from the weak lensing maps 
of D ahle et al.l (j2002l ). where, despite showing a fairly 
regular projected galaxy number density distribution, 
the weak lensing maps showed a highly elliptical mor- 
phology with major axis pointing along an axis going 
from south-west to north-eas t (similar to the X-ray 
major axis angle). iSmail et al.l ([1995) found the position 
angle for the mass, gas and projected galaxy density 
distributions within ~ 1.25 arcmin ('^ 300 kpc) are all 
aligned along the same axis. The combined ra di o and 
X-ray observations of IMazzotta fc Giacintuccil (120081 ) 
have shown a spatial correlation between the spiral 
structure and radio emission associated with a core mini 
radio halo, suggesting a population of relic electrons 
injected by the central AGN have been reaccelerated by 
turbulence associated with the gas motions causing the 
cold fronts. 

4.2.3. Abell 2142 

The first cold fronts were discovered in Chandr a 
observations of Abell 2142 (|Markevitch et all l2000f) . 
The cluster is o ptically rich (Abell richn ess class 2 ; 
Abell et all Il989l). hot (kT=9keV; Tw hitc eFall liggl 
Allen fc FabianI 119981 : iMarkevitch et aLiil998i r i2000i rX- 
ray luminous (^xfO.l - 2.4keV] = 11.4 x lO'^Vg 8"^ 
iBohringer et al.|[2000f ) with redshift z=0.089 and velocity 
dispersion a = 1280 km s'^ (|Oegerle et al.lll995l ). The 
top left panel of Figure [TOj shows the Chandra X-ray 
image of Abell 2142 which has an elliptical morphol- 
ogy with its position angle aligned along a south-east 
to north-west axis. Along this axis and to the north- 
west is the cold front for which density and temperature 
jumps were measure d, whilst a second less prominent 
cold front, noted by IMarkevitch et"al] (|2000D . lies just 
south of the core. The north-west cold front appears to 
have a cometary morphology with the sharp brightness 
edge to the north-west forming the head, and a tail point- 
ing towards the south-east where the surface brightness 
drops less sharply with radius. The temperature map 
presented in the top right panel of Figure [TU| is gener- 
ated using a much longer expos ure (Obs ID 5005 data) 
than the one shown in Markevit ch et al.l (2000) and re- 
veals a complex temperature distribution, where the two 
cold fronts clearly separate regions of cool and hot gas. 
The temperature morp hology of the central regi on is very 
similar to that seen bv IMarkevitch et al.l (|2000r ). i.e., the 
brightness peak is much cooler 6 keV) than the gas at 
larger radii, and there is a finger of ~ 7keV gas extend- 
ing from the core towards the north-west. The regions 
inside the north-west cold front are hotter than the core 
gas, at ~ 8 — 9 keV, but significantly cooler than the gas 
at larger radii outside the front (> 12keV). Note that 
towards the south-east, the gas temperature is around 
7 — 8 keV, inconsistent with the rise in temperature seen 
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in the same region bv iMarkevitch et all (|2000f ). This 
differe nce is hkely attributed to the different data sets 
used— IMarkevitch et all (|2OOO0 used early ACIS-S data 
where the cahbration was poorly known, whereas here 
the ACIS-I data are used along with the latest and much 
improved calibration files. 

The exact nature o f the cold fronts is uncertain. 
IMarkevitch et alJ ()200Clf ) propose they are due to the sur- 
vival of dense subcluster cores during a major merger, 
but later IMarkevitch fc VikhlininI (|2007[ ) claim the cold 
fronts are more likely to be the result of core gas 
motion (not unlike RXJ1720. 1+2638, MS1455.0+2232 



and Abell 1201 of this sample). However, AbeU 2142 
seems anomalous since, while there is a finger of cool 
gas extending from the core towards the north-west, 
this finger does not form a coherent spiral-like temper- 
ature structure joining the core to the north-western 
cold front, as seen in the clusters MS1455. 0+2232 and 
RXJ1720.1+2638. This of course could be due simply to 
different viewing angles, since the existence of multiple 
fronts lying on different sides of the cluster core is fairly 
compelling evidence for sloshing type cold fronts. 

Power ratio analyses of ROSAT X-ray images gave an 
indication that Abell 2142 was a dynamically active sys- 



Cold Front Clusters from the Chandra Archive 



19 




2D.Q 15.E1 10.a 

(liZDDD) 



Q5.0 14:57:Q0.D 




20.a is.a 





Fig. 9.— Same as Figure[2]but for MS1455. 0+2232 and with an SDSS r-band image at the lower left. 



tern (jBuote fc: Tsai|[l996l ). and th e case for a merger wa s 
strengthened by the analysis of iHenrv fc Briel (|1996D . 
which showed strong variations in the azimuthal temper- 
ature distribution. In an earher s tudy using 103 cluster 
member redshifts, lOegerle et al.l ()1995l ) found evidence 
for substructure at the 90% confidence level using the 
Dressler-Shectmann A test, although they were unable 
to isolate the substructure spatially or in velocity space 
and conclude Abell 2142 is a cluster with uncertain dy- 
namics. Indirect evidence for dynamical activity came 
from the observation of two dominant bright elliptical 
galaxies in the cluster core, with one having a redshift 



consistent with the cluster mean redshift and the second 
offset in peculiar velocity by ^ 1600 km s~^. Further evi- 
dence fo r dynamical activit y came from the weak lensing 
maps of lOkabe fc Umet"sul (120081 ). who found evidence 
for a significant excess in their weak lensing signal at 
a projected distance of ^ 540 kpc ahead of the north- 
west cold front. At radio wavelengths, Abell 2142 ex- 
hibits diffuse radio emission extending ~ 350kpc about 
the cluster core, although it is unclear whether this is 
part of a massive radio halo associate d with a clus - 
ter merger (Giovannini ct al. 1999). Blito n et all (|1998l ) 
show a NAT source which has a tail extending away from 
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Fig. 10. — Same as Figure[2]but for Abell 2142 and with an SDSS r-band image at the lower left. 



the cluster center and in roughly the same direction as 
the major axis of the X-ray emission. This NAT galaxy is 
displaced by ~ 1400 km s^^ from the cluster rest frame, 
and its head is located just within the edge of the north- 
west cold front (coincident with the X-ray point source 
seen to lie just within the north-west cold front in Fig- 
ure (TUl). Despite the uncertain nature of the merger ge- 
ometry, it is clear that Abell 2142 is a dynamically active 
cluster. 



5. DISCUSSION 



Selecting clusters based purely on the existence of a 
cold front has produced a sample with two distinct X- 
ray morphological types: those with clearly disturbed 
X-ray morphologies (Section 14. 1|) and those with fairly 
relaxed X-ray appearances, aside from the cold fronts 
(Section l4.2|) . In Table Hj the evidence indicating merger 
activity for the cold front sample, along with the X-ray 
morphologies, is summarized. Once again, to emphasize 
the dichotomy, the clusters in Table |4] are separated into 
the two subsamples presented in Section [4l 

Of the six clusters in Section [4T] five have significant 
evidence for merger-related activity at other wavelengths 
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and, as was evident from their X-ray morphologies, it is 
clear the existence of a cold front is good evidence for 
merger activity®. The sixth cluster, Abell 2069, does 
have e vidence for s ubstr ucture in its projected galaxy 
densitv lCioia et all (1982) , although it is not as well stud- 
ied at multiple wavelengths when compared to the other 
five clusters, thus the evidence for merger-related activity 
at other wavelengths is less compelling. Within this sub- 
sample there are a broad range of X-ray morphologies. 
This can be traced to the different formation mechanisms 
for cold fronts which depend on the details of the merger 
event (e.g. merger stage, mass of sub cluster, impact pa - 
rameter or the existence of cool cores; IPoole et al.l[2006l) . 
For example, upon comparing Abell 1201 to 1ES0657- 
558, it can be seen that both contain remnant cores ob- 
served after pericentric passage. However, Abell 1201's 
cold front is due to gas sloshing of the primary's cool 
core, while the cold front in 1ES0657-558 occurs in the 
remnant core and there is no evidence for a cold front due 
to motion of its primary's core. Two scenarios relating 
the details of the initial conditions of the mergers to the 
cold front formation can explain this: (1) The primary in 
1ES0657-558 did not contain a significant cool core, thus 
no cold fronts analogous to those in Abell 1201 are seen. 
(2) The impact parameters of the mergers differ signif- 
icantly, i.e., the merger in the 1ES0657-558 was almost 
directly head on, thus the primary's core was destroyed 
by the secondary's passage, while the pericentric passage 
of the secondary in Abell 1201 was off- axis and did not 
significantly disrupt the cool core allowing the formation 
of the cold fronts. Thus, aside from being good merger 
signposts, cold fronts also provide additional information 
which can be used in tandem with observations at other 
wavelengths to constrain the merger histories in these 

clusters. 

Turning to the clusters in Section 14.21 Table |4] shows 
that currently there is very little observational evidence 
in the literature to determine a clear cut explanation 
for the cold fronts in clusters with overall relaxed mor- 
pholog ies. The simulations of lAscasibar fc MarkevitchI 
mm) showed that pure gravitational disturbances in 
the form of infalling gasless dark matter subclumps are 
able to produce these cold fronts whilst maintaining the 
overall relaxed appearance. Commonly seen in simula- 
tions of this scenario are spiral-like temperature struc- 
tures, which we also see in the temperature maps of 
RXJ1720. 1+2638 and MS1455.0-f 2232 (Figures [8] and [H 
respectively), showing the cool gas causing the front to 
be connected to the cool cluster core. Also apparent 
in these simulations are the existence of multiple cold 
fronts at different radii; a signature clearly observed in 
Abell 2142. Marginal evide nce is fomid favoring this sce- 
nario in RXJ1720. 1-1-2638 (jOwersI 120081. where there is 
substructure in the galaxy density distribution but no 
apparent secondary structure in the X-ray image. With 
further optical data, RXJ1720.1-F2638, AbeU 2142 and 
MS1455. 0+2232 c ould provide definitive tests of t he sce- 
nario put forth bv lAscasibar fc MarkevitchI (|2006[ ) where 
an essentially gas-free group perturbs the core ICM, set- 
ting off sloshing whilst providing minimal further evi- 
dence of a merger at X-ray wavelengths. 

* Note that the existence of additional structure at other wave- 
lengths was not considered in the selection of the cold front sample. 



While much of the evidence favoring a merger origin 
for the relaxed appearing sloshing type cold front clus- 
ters comes from simulations, two significant observations 
exist which show the plausibility of merger-induced slosh- 
ing as a mechanism for the formati on of cold fronts. The 
XMM observations of Abell 1 644^ (jReiprich et al.ll2004t 
iMarkevitch fc Vikhliiiiiil 120071) and the Ch andra obser- 
vations of Abell 1201 presented here and in lOwers et al.l 
(2009b|) both reveal sloshing cold fronts and also rem- 
nant cores which are apparently causing the perturba- 
tion. The optical follow up observations of Abell 1201 
showed evidence for localized yelocity substructure coin- 
cident with the remnant core (jOwers "eraI1l2009bn . con- 
firming it is a significant merging substructure, although 
the latest optical observations of Abell 1644 do not re- 
veal significant evidence for localized velocity substruc- 
ture (jTustin et al.l 12001*). Without detailed simulations 
of these two particular objects, it is unclear whether 
they will evolve into relaxed appearing clusters with cold 
fronts, like those presented in Section l472l Nevertheless, 
it is clear that both Abell 1644 and Abell 1201 provide 
strong empirical evidence of merger-induced sloshing as 
a viable mechanism for producing cold fronts. 

Perhaps the most appealing outcome that studies of 
cold front clusters offer is the possibility of using them 
to distinguish major and minor mergers. Clearly the 
clusters in Section 14.11 are systems which are undergo- 
ing significant dynamical evolution and are most likely 
major mergers. However, the same conclusion can- 
not be drawn based on the evidence available for the 
clusters presented in Se c tion 14.21 As discussed above, 
lAscasibar fc MarkevitchI ()2006l ) found that only a rela- 
tively minor merger was required to produce cold fronts 
in relaxed appearing clusters. Future dynamical analyses 
of RXJ1720. 1+2638, AbeU 2142 and MS1455.0+2232 at 
optical wavelengths will determine the level of dynami- 
cal activity required to produce a cold front. If it turns 
out to be true that the appearance of these clusters can 
be explained by minor mergers, then the combination of 
the existence of a cold front and X-ray morphology will 
provide a useful diagnostic for discriminating major and 
minor mergers in clusters. 

6. SUMMARY AND CONCLUSIONS 

We have presented and described in detail a sample 
of nine clusters selected from the Chandra archive based 
purely on the existence of cold fronts in their ICM. These 
clusters are the subject of an ongoing study to deter- 
mine if cold fronts may be used to reliably identify recent 
cluster-cluster mergers. 

In selecting the cold front sample, particular atten- 
tion was paid to obtaining highly reliable identifications. 
To this end, each cluster met the following strict selec- 
tion criteria: A total Chandra exposure time greater than 
40ks, 0.05 < z < 0.3, a visible surface brightness edge in 
the X-ray image which has a density jump exceeding 1.5 
at the lower 90% confidence level, the measured gas tem- 
perature must be lower on the denser side of the front, 
and the pressure must be continuous across the front 
within the errors. This resulted in the identification of 
nine cold front clusters for which we present temperature 

^ AbeU 1644 is not part of the cold front sample presented in 
this paper, since its redshift falls below the cutoff of z=0.05. 
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maps, X-ray and optical images and an overview of the 
multiwavelength evidence for merger activity available in 
the literature. The key conclusions of this study are as 
follows: 

• the measurement of thermal properties across the 
fronts was an essential diagnostic tool; Abell 2034 
and Abell 665 were rejected on this basis, since 
their fronts were better explained as shocks. These 
are very likely merger shocks and if these clusters 
were to remain in our sample, the contamination 
would not be harmful when considering the aim is 
to use cold fronts as reliable signposts of merger 
activity. 

• There is a dichotomy within the sample of cold 
front clusters: those which harbor clearly disturbed 
X-ray morphologies, and those which appear to 
have fairly relaxed X-ray morphologies, apart from 
the existence of a cold front. 

• The majority of those clusters with disturbed X-ray 
morphologies have significant evidence for merger 
activity at multiple wavelengths, while minimal ev- 
idence for a merger exists at other wavelengths in 
the relaxed appearing clusters. 

• Two of the relaxed appearing clusters 
(RXJ1720.H-2638 and MS1455. 0-1-2232) ex- 
hibit spiral-like structures in their temperature 
maps, similar to those seen in the simulations of 
lAscasibar fc MarkevitchI (|2006D . 

To conclude, in the case of the cold front clusters with 
disturbed morphologies, it is clear that a cold front is 
an excellent signpost for ongoing major merger activity. 
However, it appears verification that all cold fronts are 



signatures of merger activity (major merger or otherwise) 
relies critically on the results of follow up observations of 
the relaxed appearing cold front clusters. The sample 
is the subject of extensive follow up observations using 
multi-object spectroscopy at optical wavelengths to ob- 
tain large samples of confirmed cluster members. These 
observations will be used to search for dynamical sub- 
structure, constrain merger dynamics and to study the 
effects of the merging environment on the cluster galax- 
ies. 
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APPENDIX 

SURFACE BRIGHTNESS BROKEN POWER LAW DENSITY MODEL 

For the purpose of quantifying the edges in surface brightness, it is assumed that a discontinuity is formed at the 
interface of gases of differing physical properties, occupying two similar concentric spheroids with their rotational 
symmetry axis in the plane of the sky. Thus the curvature of the front relative to the line of sight is assumed to be 
the same as that in the plane of the sky. The electron density, ng, is described by the broken power law function 



rieir) 



ne,i 

ne,2 



it)' 



r<Rf, 
r>Rf, 



(Al) 



where Rf is the radius of the discontinuity, ne,i and rig, 2 are the densities of the inner and outer gas at i?/, respectively, 
the elliptical radius, r, is defined by = -I- e^^C^, C is the coordinate along our line-of-sight, t,^ determines the axis 
length in that direction and ro, the elliptical radius in the plane of the sky, is defined by 



w"^ = {x cos 9 + y sin ( 



''{ycosi 



X sm ( 



(A2) 



where x and y are cartesian coordinates with respect to the center of the ellipse, is the position angle and e is 
the ellipticity. The X-ray emission primarily arises from Bremsstrahlung radiation with secondary contributions from 
collisional line excitation radiation. Both these emission processes depend on the electron and ion densities in the 
ICM, and the X-ray surface brightness, Sx, is proportional to the integrated square of the electron density along the 
line of sight so that 

^^^mttif/^^^'^^"^''^^' ^^^^ 

where z is the cluster redshift, rig is the electron density, and F{T, Z) gives the dependence of the telescope response to 
thermal emission from an optically thin gas of the temperature, T, and abundance, Z, (there is also some dependence 
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on the column density, n/f , and z) integrated over the energy range of interest (0.5-7 keV). The function F(T, Z) also 
contains the constant of proportionality relating Ue to the proton density. Substituting the relations from Eqns. lAll 
into Eqn. IA3I and integrating, the surface brightness distribution is 

where (3j^i — (3{^,ai — 5; 1 — (j^)^)' f^i-^ ~ f^ih^^"^ ~ l'^^ ~ ("^)^) incomplete beta functions, and P2 = 
P{^,ce2 — ^) is a complete beta function. For ai < 0.5, recursion relations were used to overcome the divergence of 
the incomplete beta function. At ai — 0.5, —0.5, the integrals were solved analytically. The constants Ai and A2 are 

_ nl,F{T,,Z,)Rf _ nl2FiT2,Z2)Rf 

- 4^(l + z)4e^ ' - 4^(l + z)4ec ' ^ ^ 

and the density jump at the discontinuity can be measured by taking the ratio of Ai to A2 



_ AiF{T2,Z2) 



(A6) 



ne,2 yA2F(Ti,Zi) 

For temperatures greater than 2keV, F{T,Z) is quite insensitive to temperature. Thus, the density jump is well 
approximated by « \[^^' 

Note that while the density is discontinuous at the front, after projection onto the sky, the surface brightness is not. 
Thus we refer to the corresponding feature as an "edge" in the surface brightness. As illustrated in Figure [1] surface 
brightness does change abruptly at cold fronts, which can give the misleading impression of a discontinuity. 
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TABLE 4 

Summary of merger, evidences and X-ray morphologies for the cold front sample. 



Cluster 


X-ray Morphology 


Temperature Map 


Optical Substructure'' 


Radio 


Lensing-based 


References 






Structure 






mass maps 






Disturbed Clusters 


1ES0657-558 


disturbed, remnant core 


shock front, 


projected galaxy density. 


halo 


substructure 


1, 3, 10 






cool remnant core 


position plus velocity 








Abell 3667 


disturbed 


cool mushroom-cap. 


projected galaxy density 


reUc, NATs 


substructure 


9, 17, 18 






hot patches 


position plus velocity 








Abell 1201 


remnant core, elliptical. 


cool finger joining 


projected galaxy density. 






16 




cool core 


cold front and cool core position plus velocity 








Abell 2069 


elliptical, remnant core 


cool remnant core. 


projected galaxy density 






6 






hot patches 










Abell 1758N 


disturbed, bimodal 


cool remnant core 


projected galaxy density 


halo, NATs 


substructure 


4, 8, 12, 14 


Abell 2163 


disturbed, remnant core 


cool remnant core. 


projected galaxy density 


halo, NATs 


substructure 


5, 11, 15, 19 






hot patches 


position plus velocity 








Relaxed Appearing Clusters 


RXJ1720.1+2638 


relaxed, cool core 


cool spiral, cool core 


projected galaxy density. 




possible substructure 


4, 15 






hot patches 


possible position plus velocity 








MS1455.0+2232 


relaxed, cool core. 


cool spiral, cool core 






elliptical distribution 


4 




mildly elliptical 












Abell 2142 


elliptical, cool core 


cool core, cool finger 


possible position plus velocity 


NATs, 


substructure 


2, 7, 13, 14 






hot patches 




possible halo 
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^ Position plus velocity substructure refers to substructure found in tests which combine spatial and velocity information. Projected galaxy density substructure refers 
to bimodality or multimodality in projected galaxy density maps. 



